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Abstract. A new formalism is introduced for the trans- 
fer of polarized radiation. Stokes parameters are shown 
to be four-vectors in a Minkowski-like space and, most 
strikin gly, the radiative transfer equation (RTE) turns out 



The fundamental problem addressed in the present 
paper is the existence of an analytical solution to the 
RTE. At present, analytical solutions exist for o nly very 
limited cases and the known formal solutions ( Landi 



to be i n infinitesimal transformation under the Poincare DegPInnocenti and Landi Degrinnocenti, 1985 ) do not 



(plus dilatations) group. A solution to the transfer equa- 
tion as a finite element of this group is proposed. 

Key words: Radiative transfer - Polarization 



1. Introduction 



Since the pioneering paper of Unno (1956) in which for 
the first time a transfer equation was derived for polar- 
ized light in the presence of a magnetic field, a num- 
ber of attempts have been made to obtain a solution for 



it. Rachkowksy (1967), after completing the equation by 
adding anomalous dispersion effects, gave the first analyt- 
ical solution for the case of Milne-Eddington atmospheres. 
While numerical solutions were more and more successful 



(Wittmann, 1974; 


Landi Degl'Innocenti, 1976; 


Rees et al.. 


1989 


; Bellot Rubio ct al., 199^; 


Lopez Ariste and Semel, 


199C 


) analytical trials struggled to overcome the problem 



of non-commuting absorption matrices. This problem, al- 
ready noted in a pap er of Landi Degl'Innocenti & Landi 
DeglTnnocenti (1985), is pointed out clearly in the work of 
Semel & Lopez Ariste ( 1999[ ) (hereafter referred as Paper 
I) as the origin of all previous limitations. In this paper, 
deep insights are given on how to benefit from the physical 
significance of the various terms in the transfer equation 
to attain an analytical solution as general as possible. Sev- 
eral transformations used here to simplify the RTE seem 
to indicate that the physics of the Stokes parameters is 
best described in geometrical terms. Indeed, we shall show 
in this paper that the Poincare (plus dilatations) group is 
at the origin of these geometrical aspects, and we shall 
use its algebraic properties to solve the problem even for 
non-commuting absorption matrices and give a general 
solution for this equation of transfer. 
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offer, in the general case, any advantage from the compu- 
tational point of view. Apart from the interest of a com- 
putable analytical solution by itself, it would be of greater 
importance for testing numerical codes which nowadays 
are considered acceptable only by observing their con- 
vergence upon increasing number of layers. This is not 
a very satisfactory situation. We think that the solution 
presented in this paper is the first step towards a general 
solution suitable for testing computations. 

Why an analytical solution could not be found so far is 
our first question. It is our belief that group theory is the 
key to the solution and that explains why it had not been 
reached up to now (this could already be understood from 
the conclusions in Paper I which stressed the importance 
of the non-commutativity of matrices in this problem, but 
it is still clearer here) . Group theory is not common in as- 
tronomical literature. However the RTE of polarized light 
has become a more and more important problem in astro- 
physics, related, for instance, to the measurement of mag- 
netic fields via the Zeeman or Hanle effects. Hence group 
theory should be accessible to the concerned astrophysical 
community, and with its help we can give a method to find 
the solution, and explicitly give the full expression of the 
analytical solution for the most general case. Any other 
particular or general solution will benefit from the use of 
advanced linear algebra, and so avoid wasted effort. Last 
but not least, it may deepen our understanding of Stokes 
polarimetry. 

We begin our research by disclosing the mathematical 
nature of the Stokes vector, starting with its physical def- 
inition and extending to the appropriate mathematics to 
treat our problem. Extensive literature has already been 
devoted to the existing relations between polarization and 
the Lorentz group, mainly from the optical point of view 
(see for instance Cloude 1986; Givens & Kostinski 1993; 
Sridhar & Simon 1994 and references therein). Usually 
these works take off from the Jones formalism for polar- 
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ized light (Jones, 1941) and develop these relations. Here 
a similar path is followed to show that the Stokes vector 
is a 4-vector in a Minkowski-like space. The demonstra- 
tion is based on the comparison between the usual defini- 
tion for the Stokes parameters (see for example Shurcliff, 
1962, or Jefferies, Lites & Skumanich, 1989) and the well- 
known relations between the definition of a spinor and its 
different representations (see for instance Landau & Lif- 
shitz, 1971). In the same line of thought, we will propose 
in Section 3 that the RTE is just a representation of an 
infinitesimal Poincare (plus dilatations) transformation in 
the Minkowski-like space where the 4-vectors are best de- 
scribed. This new interpretation of the equation of transfer 
suggests that any solution to this equation must be a fi- 
nite Poincare transformation. We calculate it in Sections 
4 and 5. We note that nothing differentiates the trans- 
fer equation as written elsewhere in the literature from 
a Poincare transformation, although at present a demon- 
stration is not available, apart from this similarity. We 
anticipate that a finite Poincare transformation may be a 
solution for the RTE. 

The concepts of 4-vector, Minkowsky space or 
Lorentz and Poincare transformations must be understood 
throughout this paper in their purest mathematical sense, 
beyond their historical meaning. These concepts origi- 
nated in the framework of the special theory of relativ- 
ity, however mathematics did abstraction of these tools 
and incorporated them into more general frames of ge- 
ometry and group theory. Therefore we define 4-vectors 
as sets of 4 numbers characterized by their Minkowsky 
norm and described in a hyperbolic 4-dimensional space 
called the Minkowsky space, which, in this paper, we will 
refer to as the Minkowsky-like space, in order to stress 
the difference with the usual Minkowsky space used in 
relativity. Lorentz and Poincare transformations describe 
movements in this space: generalized translations and ro- 
tations. In the relativistic formalism these movements are 
interpreted as a change of reference system. In this paper 
they are not given that meaning, but are seen as changes 
in the polarization state. The manipulation of these con- 
cepts is identical here and in the relativistic formalism, 
if one takes care in substituting the four Stokes parame- 
ters for space and time, and for the speed its analogues as 
shown in Section 3. 



2. Stokes parameters as a 4 vector 

The 4 Stokes parameters are usually represented by 
/, Q, U, y, where / stands for the total intensity of light, Q 
and U for the linear polarized light in two axes rotated by 
45° one from the other, and V for the circularly polarized 
light. These four quantities are not completely free but 
must satisfy an energy condition: there cannot be more 
polarized light than total light. This condition suggests 
an interpretation of the Stokes parameters as a 4-vector 



in a Minkowski-like space. The norm of a vector in such a 
space is defined as 



(1) 



Hence the above condition is naturally satisfied by vectors 
with a positive norm, in parallelism to time-like vectors in 
special relativity. Continuing this parallelism, a light-cone 
can be defined as the surface which obeys the condition 



11^11=0, 
that is 



(2) 



(3) 



In our context, this surface contains all the different pos- 
sibilities for fully polarized light. The Stokes vector must 
be inside or, in the limit, on this light-cone (see Fig. ^ to 
obey condition (^. An exception to this parallelism with 
special relativity: the backward light cone does not have 
an equivalent with the Stokes vectors. 




Fig. 1. 3-dimensional projection of the light cone for the 
Stokes vector 



To consolidate and extend this interpretation of the 
Stokes parameters we begin with the definition of the 
Stokes parameters in terms of the electric field. 

A transversal monochromatic light wave is completely 
described by and Ey, the components of the electri- 
cal field in a plane perpendicular to the direction of the 
propagation of light, z. Following definitions in Landau & 



Lifshitz (1971), these two components can be arranged in 
a 2-dimensional vector. Since it transforms linearly under 
the proper Lorentz group (ibid.) this vector can be called 
a spinor of rank one. As an illustration of this kind of 
transformation, we profit from the fact that any element 
of this group in its 2-dimcnsional representation can be 
written as a linear combination of the Pauli matrices (plus 
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the 2x2 identity matrix): 



(^2 



1 

1 

1 

1 



CTi 



0-3 = 



1 

-1 

-i 

1 



(4) 
(5) 



and transform the electric field vector by these matrices: 



0-2 



CT3 



Ex 

Ey 

Ex 



Ex 

-Ey 
-iEy 

iE^ 



(6) 



The result is always a new 2-dimensional electric field vec- 
tor which describes a different state of polarization. 

A spinor of rank two can be easily constructed by mul- 
tiplying conveniently two spinors of rank one. For our par- 
ticular spinor we obtain 



J' 



TP Tr* TP TP* 



EyE* 



EyEy 



where the * symbol stands for complex conjugated. This 
2nd rank spinor is to be compared with the coherency 
matrix (see Born & Wolf 1980). In fact the last is defined 
for any given light beam, and a mean over frequencies or 
time is necessary in the above expression of J' to fulfill the 
definition. The average of J' is a linear combination of ma- 
trices of the form referred. Since spinors of rank two form a 
linear vector space, any linear combination of spinors will 
be a spinor as well. Hence the coherency matrix, defined 
usually as 



< ExE* 

< EyE* 



> 
> 



< ExE; > 

< EyE* > 



is still a spinor of rank two. 

Since a spinor built in this way has 4 independent com- 
ponents (the four entries of the matrix), there must exist a 
relation between it and a 4-vector, which also has 4 inde- 
pendent components. Technically speaking, both must be 
different realizations of the same irreducible representation 
of the Lorentz group (Landau & Lifshitz 1971, page 55). 
The components (/, Q,U,V) of this 4-vector are indeed 
related to the components of J as 

1 



Q 



U 



V = 



^(Jr2 
1 



J21), 

-J21). (7) 
These components are identical to the definition of the 



Stokes parame ters a s given by Jefferies et. al. (198£), or 
Born & Wolf (|l980|) for instance. The conclusion is evi- 
dent: the resulting 4-vector, derived from the coherency 



matrix, a spinor of rank 2, is the Stokes vector. This for- 
malism can be given in an alternative way: an usual basis 
for spinors is the set of Pauli matrices plus the 2x2 iden- 
tity matrix. The coherency matrix expressed in this basis 
has for coefficients the Stokes parameters 



J = /(To + Q<^i + Ua2 + Vas 



(8) 



Using this relation we can write in a more compact form 
relations (0) as: 



lTr[^a], 



(9) 



where the vector a has four components, the 2x2 identity 
matrix and the 3 Pauli matrices, and Tr denotes the trace 
operation on matrices. 

These relations stress further the interpretation of the 
Stokes parameters as a 4-vector in a Minkowski-like space. 
We stress that when talking about a Minkowski- Zifce space 
we mean that the coordinates in our 4-dimensional space 
are no longer space and time, but the Stokes parameters, 
contrary to the usual Minkowski space used in relativistic 
formalism. On the other hand the underlying hyperbolic 
geometry is exactly the same in both cases, mathemati- 
cally speaking they are the same space. All the usual prop- 
erties of the Stokes parameters are recovered in this space. 
For instance the sum of two 4-vectors is a new 4-vector, a 
well-known property of the Stokes parameters. Contrary 
to the usual Minkowsky space, the absence of the back- 
ward light cone implies not only that negative intensities 
are meaningless, but also that the negative of a Stokes 
vector does not exist. Hence subtraction of Stokes vectors 
is naturally forbidden in this space. 

3. Radiation transfer equation as a Poincare (plus 
dilatations) infinitesimal transformation 

Since the Stokes vector is a 4-vector in a Minkowski-like 
space, one may wonder what would be the meaning of a 
Lorentz transformation over the Stokes parameters. 

Homogeneous Lorentz transformations form a 6- 
parameter Lie group: 6 generators suffice to describe all 
possible infinitesimal transformations. These generators 
are (see for example Greiner 1990, or any textbook in spe- 
cial relativity or group theory)]^ 

— The three 4x4 matrices for 3-dimensional spatial ro- 
tations 
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^ Everywhere in this paper we use g — diag{l, —1, —1, —1) 
as the metric for the Minkowski space 
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The three 4x4 matrices for hyperbolic rotations (or 
Lorentz boosts in relativistic terms) 



Kn = 



/O 1 0^ 

10 



Vo 0. 



/O 1 0^ 



10 

Vo 0. 



^0 r 





Vi 0. 



An infinitesimal Lorentz transformation over the Stokes 
4-vector / can be expressed as a sum of generators mul- 
tipHed by their respective infinitesimal parameters: 



1=1 



E 

i=Q,U,V 



(10) 



i=Q,U,V 



For convenience, we can re-write all these 
parameters(/?i, 7i) in terms of a common parameter 
dr, expressed in differential form (we are dealing with an 
infinitesimal transformation) 



/3 = (/3Q,/3[/,/3y) = -{pQ,pu,pv) ■ dr, 



7 = ( hq.yrT.^v) = -ivq.Vjj.Vv) ■ dr. 



(12) 



And the infinitesimal transformation reads now: 

\i=Q,Uy i=Q,U,V J 

Given the infinitesimal character of the transformation, 
expressed explicitly by putting the common parameter dr, 
it easily leads to a differential equation for / in the variable 
r: 



— I = -K'I, 
dr 

where K' is a 4 x 4 matrix given by 

K'= p,s,+ 

i=Q,U,V i=Q,U,V 



(14) 



/ 

Vu 
Vv 







w 





Pv 


-pu 


-Pv 





PQ 


Pu 


-PQ 






Matrix K' looks like the well-known absorption matrix 



(Landi Degi'Innocenti, 1992). At this point the possibility 



that the RTE could be written as an infinitesimal transfor- 
mation involving the Lorentz group seems to be at hand. 
But K' still differs from the general form of the absorption 
matrix. In particular the inclusion of a diagonal term is 
necessary if we want to take into account the scalar ab- 
sorption represented by 77/ in the usual matrix. To include 
it we just add to the usual 6 generators of the Lorentz 



group the one for the dilatation transformation ^. While 
dealing with the homogeneous group, we can represent the 
generator of dilatations by the 4x4 identity matrix. We 
will repeat all the previous steps calling the new parame- 
ter for this transformation ?// : 

\ i=Q,U,V i=Q,U,V ) 



to finally obtain: 
where K is given by 



(15) 



+ 


E 


PiSi + 


E 


i- 


=Q,uy 


i- 


=Q:Uy 


VI 


VQ 


Vu 


Vv 


VQ 


VI 


PV 


-pu 


Vu 


-pv 


VI 


PQ 


Vv 


Pu 


-PQ 


VI 



Matrix K can now be compared to the so-called absorp- 



(11) tion matrix which appears in RTE (Landi Degi'Innocenti 



1992; Jefferics et al., 1989). Coefficient rji gives the scalar 
absorption, independent of polarization. This absorption 
has its equivalent in a contraction (that is a negative di- 
latation) of the Stokes 4--vector. The 3-vector rj is respon- 
sible for the creation and absorption of polarization, which 
is understood here to be a hyperbolic rotation (or Lorentz 
boost) of the Stokes 4-vector. Finally the 3-vector p gives 
the so-called Faraday rotation in Zeeman effect and, as 
its name seems to indicate, it rotates the Stokes 4-vector 
inside the 3-dimensional space of polarized states (also 
called Poincare sphere). Note that from the previous para- 
graph it cannot be stated that any infinitesimal Lorentz 
(plus dilatations) transformation is a transfer equation. 
We state that the absorption matrix as it is usually written 
elsewhere in the literature, cannot be differentiated from 
an infinitesimal Lorentz transformation (indeed we have 
decomposed the absorption matrix in the RTE in terms 
of the infinitesimal generators of the Lorentz transforma- 
tions). In fact, it is well known that the actual absorption 
matrix for Zeeman effect must still obey some further con- 
straints. In this sense a general Lorentz (plus dilatations) 
transformation is too general: absorption matrices stand- 
ing for a real physical process would constitute only a 
subset of all possible Lorentz (plus dilatations) transfor- 
mations. While it is evident that a further study of this 
relation is necessary, in this paper we will only make use 



^ The dilatation transformation does not belong to the usual 
Lorentz group. A usual definition of the transformations be- 
longing to this group is that they do not change the Lorentz 
norm of any 4-vector. By definition, a dilatation transforma- 
tion does change this norm. Fortunately the new set of 7 gen- 
erators is still a group. 
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of the mathematical advantage and postpone the rest for 
a forthcoming paper. 

Eq. ( p^ is not yet the complete transfer equation. An 
inhomogeneous term, the emission vector, is still needed. 
For this purpose there are appropriate movements in the 
Minkowski space: the homogeneous Lorentz group can be 
extended to the inhomogeneous Poincare group. This 10- 
parameter Lie group shares 6 infinitesimal generators with 
the Lorentz group and adds 4 more generators (P/, Pq, Pu 
and Py) to take into account translations along I, Q, U 
and V (what in relativistic formalism would be transla- 
tions in time and space). These generators operate on a 
generic Stokes vector as follows: 

/1\ /o- 

(16) 



Prl 



Pul 







.Pvl 



After including dilatations, the infinitesimal inhomo- 
geneous transformation is given by 



1=1- KIdT + dr 



=I,Q,U,V 



I - KIdT + JdT{n) 



where J is the emission vector. Some algebraic manipula- 
tions equivalent to those used for the homogeneous group 
will lead to the complete radiative transfer equation. 

It is interesting to note the fact that these generators 
allow one to write an inhomogeneous term (the source 
function) in a pseudo-homogeneous way: 

\ i=I,Q,U,V ) 

A useful representation to understand this apparent para- 
dox is the one using differential operators. Let us call 
// = /, Iq = Q and so on, and write 

where i = I,Q, U, V. It is evident that with such an oper- 
ator, relations ( |l6| ) hold. All the other generators can be 
rewritten in this representation. For instance, the dilata- 
tion generator can be written 

d 



D 



E 



i=I,Q,U,V 

and the generator of rotations in the plane UV, for in- 
stance, is 



So = -Ii 



d 



d 



div diu ' 
The interested reader will find good discussions on the 
representations of the Poincare group and the dilatation 
transformation in Greiner (1990), Gourdin (1982) or Jones 
(1996) for example. 



4. Finite transformations as a solution 

In view of the results obtained in the last section, the 
solution to the RTE appears to be quite straightforwardly 
a finite Poincare (plus dilatations) transformation. The 
important fact now is that we already know how to write 
such a finite transformation: if we denote the 11 generators 
by ti, a finite element of this group can always be written 
as 



exp 



where £,i are the parameters of the transformation for each 
movement, the equivalent of angles for usual rotations. 
The next problem is how to calculate those finite parame- 
ters from their infinitesimal counterparts ( ?]/, rjq and so 
on). Unfortunately this is not an easy task. The problem 
resides in the non-commutativity of the generators. The 
next paragraph proposes a solution to this problem, but 
first we consider it usef ul to clarify why other approaches 
will not work. Magnus ( 1954 ) (a brief introduction to this 
paper can be found in the Appendix of the paper by Semel 
& Lopez Ariste, 1999) has already written a finite trans- 
formation in terms of such a unique exponential. A quick 
inspection of the expression given there illustrates why we 
consider that this calculation is not an easy task. Nev- 
ertheless some simple cases can be proposed in which the 
relation between the and the respective infinitesimal pa- 
rameters is plain. For instance in the case of 3-dimensional 
rotations it is always possible to transform (by means of 
the Euler angles) our initial reference system into another 
one for which the rotation axis is parallel to one of the 
new reference axes. In the case of a fixed axis, just one 
generator suffices to describe the movement. If we call this 
generator t, and the infinitesimal parameter c?^, the finite 
transformation results in 



exp (^t) = exp ^t j d£_ 



Apparently, our problem is solved if we deal with a unique 
generator per exponential. Hence, we can propose for in- 
stance a solution in the form of a product of several expo- 
nentials, one for each generator: 

Y[ exp(^,t,). 

1=1,11 

Now the derivative of each exponential in the product can 
be easily calculated in a very compact form, and so the 
product of all of them, and consequently the new finite pa- 
rameters S,i (as in the above example). But a new question 
arises: which order should be chosen for the exponentials? 
Again due to the non-commutativity of the infinitesimal 
generators, different orders produce different results, 



(18) 



and in general no particular order will be the solution. 
An answer to this kind of problem has been given by Wei 
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& Norman (1963). In what follows, we sketch the solution 
there proposed and apply it to our particular problem. We 
start with the homogeneous equation (just the Lorentz 
group plus dilatations) and in the next section we will 
incorporate the inhomogeneous part and handle the full 
Poincare group. 

The ordering problem can be recast as follows: we may 
say that by introducing an order in the exponentials we 
introduce an error. In spite of this error, let us choose a 
particular order for the exponentials and then substitute 
the (which could be calculated straightforwardly as the 
integrals over the path of the corresponding infinitesimal 
parameters) for some unspecified scalar functions gf. 



Y[ exp (giU) 



(19) 



i=l,7 



The new functions gi have to, in a certain sense, take into 
account the effect of the and correct the error introduced 
by the chosen order. Existence for those gi functions can 
only be ensured after introducing the proposed solution 
into the RTE. The following consistency equation for the 
gi's is obtained as a necessary condition for Eq.(p^) to be 
a solution: 



1=1,7 



JJexp(gjtj) 

i=i 



j=i-l 



where the dot denotes derivative over the integration vari- 
able T. This equation was the aim of this section so far. A 
less heuristic but more direct way of introducing it is to 
look for a solution of the form ( p^ and introduce it into 
the transfer equation. It is straightforward to fall upon 
Eq.(pO|) as the condition for ( p^ ) to be a solution. 

In what follows in this and the next sections we shall 
solve Eq.(|20|) for the ^^'s. The details being quite technical, 
the reader may wish to skip this and go directly to Eq. 

Eq. ( |20| ) is quite an involved equation. It requires cal- 
culation of a non negligible number of expressions of the 
form 

exp(ffitj)t, exp{-gjtj). 

This is to be done by means of the Baker-Hausdorff for- 
mula, which states that 



Y+[X,Y] + ^[X, [X,Y]] + ^[X, [X, [X,Y]]]- 



[Si, Sj] 
[K., 11] 



EyfeSfe 


0, 



which does not facilitate calculations at all. A suitable 
combination of generators will yield a new base of gener- 
ators with a gentler (from our point of view) Lie algebra. 
For instance, the following set. 



Hi 

H2 

H3 

L,; 



(Ky-Sv)-i(Ky + S[/), 

(Ky + Sv) + i(Kv -S[/), 
H* for i = 1,2,3, 



to which we add the identity, 1, for completion, obeys the 
following commutation rules 



[Hi,H2] 


= 2zH2, 


[Hi,H3] 


= -2iH3 


[H2,H3] 


= -4iHi 


[Li,L2] 


= -2iL2, 


[Li,L3] 


- 2iL3, 


[L2, L3] 


= 4zLi, 




= 0,Vz,j 



The initial algebra of 6 generators has been decomposed 
into two sub-algebras of 3 generators each, with the partic- 
ularity that each generator of one sub-algebra commutes 
(mh every generator of the other one. The dilatations gen- 
erator, which must be added to them, commutes with ev- 
ery other generator (remember that, while constrained to 
the homogeneous group, the dilatation generator can be 
represented by the identity matrix) and therefore there is 
no ordering problem associated with it. The initial prob- 
lem of ordering 7 exponentials is reduced to ordering a 
subset of 3 exponentials; the order of each one of the H's 
with respect to the L's generators or the identity being 
immaterial. 

We can rewrite solution (|l^) using the new set of in- 
finitesimal generators in explicit form as 



I[t) 



gS3(T)H3gg2(-r)H2ggi(T)Higg6(-r)L3gS5(T)L2g34(-r)Ligg7(T)]lj^^^-j^ (^21^ 



In the form they have been written in the last section, 
the generators of the Lorentz group obey the following Lie 
algebra: 



where a very special order has already been chosen. Difi'er- 
ent orders, while yielding equivalent solutions, can make 
calculations affordable or desperate. With this problem in 
mind we have chosen a particular order. 
We now express K in the new basis: 

K = ^a,H,; +^6,U 77711 

i i 

where 



ai = -^{PQ + im) 

0-2 = ^[(W + Pv) + i{riv 



Pu)] 



as = -^[{ilu - pv) - i{riv + Pu) 



b, = 



,Vi = 1,2,3. 
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Calculation of the Baker-HausdorfF scries for the new gen- 
erators simplifies a lot, and Eq. ( pO| ) reads now: 

- ^ fljHj - ^ biU - = grl + 53H3 + 

i i 

+92 [H2 + 4i.g3Hi ~4g2H3] + 

+gi [-2i52H2 + (1 + 85253) Hi + (2^53 + %iglg2) H3] + 

+561-3 + 55 [L2 - 4z56l-i - 45IL3] + 

+54 [2Z55L2 + (1 + 85556) Li - (2^56 + 815^55) L3I22) 

Comparing coefficients on both sides, and after proper re- 
arrangement the following set of differential equations is 
obtained: 

51 = -ai -I- 4z53a2 

52 = -2*5201 - (f -I- 85253)a2 

53 = 2z53ai -I- 45^02 - 03 

54 = -bi - 4156^2 

55 = 215561 - (1 + 85556)62 

56 = -2igebi + 45562 - 63 
97 = -VI 



All of them share the same boundary condition, namely 
9i{'To) = 0, to satisfy the boundary condition of the RTE. 
The function 57 can be integrated at once to give 



57(t) 



Viir')dT'. 



As expected, because of the particularities of the Lie al- 
gebra, the set of equations for 51,2,3 is separated from the 
one for 54,5,6, and each set is the complex conjugated of 
the other, so that the solution to 54.5,6 is straightforward 
once the one for 51,2,3 is given. Furthermore each set can 
be solved by quadrature, as equations for 51 and 52 de- 
pend only on 53 , whose equation is disentangled from the 
others: 



53 = "53 + P93 + li 



(23) 



where we have defined a = 4a2, l3 — 2iai and 7 — —03. 
This is a Riccati equation, and for its solution the explicit 
dependences of 01,2,3 on the integration variable r are re- 
quired. For a constant K matrix the solution is straight- 
forward, and from it those of 52 and 51. More complex 
dependences must be carefully managed (see for example 
Carihena & Ramos (1998) and references therein for the 
integrability conditions of the Riccati equation). 

Before passing to the next section, where we will gener- 
alize the method to the full Poincare group, we go back to 
Eq.(^l|). Once we have integrated the Riccati equation and 
obtained all the 5i's, we still need to calculate the expo- 
nentials. To this end we profit from a remarkable property 
of matrices H, and L^: 

Hl = Ll = -1 
Hl = Hl = Ll = Ll = 0, 



by means of which: 

e^i^i = cos5i]l + sin5iHi 



gS4Ll 

395 L2 



cos 54! + sin54Li 

]1 + .95L2 
]1 + .96L3. 



The final complete solution for the homogeneous part re- 
sults in 

lir) = []l + 53(T)H3]-[]l + 52(r)H2] 
[cos5i(T)]l-fsin5i(T)Hi] 
•[]l + 56(r)L3]-[]l + 55(T)L2]- 
[cos54(t)]1 + sin54(T)Li] • exp[57(T)]l]7(ro). (24) 

The validity of this solution is almost evident: Its deriva- 
tive results in the RTE just by making use of the differ- 
ential equations satisfied by the functions 5^. 

This is a solution to the homogeneous equation 

dr 

Let us write this solution as 

Il{t) = 0(t,to)Il(to), 

where the explicit form of 0(r, tq) can be found by com- 
paring this expression with the complete one in Eq. (p^ . 
This operator 0(r, tq) is often referred to as the evolu- 
tion operator (see mainly Landi DeglTnnocenti & Landi 
DeglTnnocenti 1985, who first introduced it). This opera- 
tor trivially obeys the homogeneous equation 



dr 



0(t,to) = -KO(T,ro), 



(25) 



with initial condition 



0(To,ro) = 1. 

Eq. ( p^ provides on its own a general analytical solution 
for the evolution operator. 

This solution is a fully general expression for finite 
Lorentz transformations plus dilatations. But radiative 
transfer cases do not cover the full spectrum of Lorentz 
transformations. In this sense the obtained solution is too 
general, in agreement with Section 3. As an illustration, 
consider the case when the p's and Ty's are zero except for 
TjQ. The function 51 becomes 



5i 



and a term of the form 



Vqdr, 



cos 5i — cosh — j r/gdr 



8 



A. Lopez Ariste & M. Semel: Analytical solution to RTE 



appears in the final solution. The hyperbolic cosine grows 
monotonously with its argument, therefore the intensity 
of the out-coming light would grow also monotonously for 
a semi-infinite atmosphere. This is a completely nonsen- 
sical result. To recover physical sense one must impose 
some constraint on the allowed transformations. This con- 
straint evidently imposes a relation between rjQ and 77/, 
whose explicit form is outside the scope of this paper, but 
which should be derived from the assumed physical pro- 
cesses. This example can be extrapolated to all the 77's and 
p's. The relations thus obtained will constrain the Lorentz 
transformations to a subset of matrices for which, never- 
theless, the above solution (|2^ ) will remain valid. 

5. Solution for the complete inhomogeneous 
equation 

To solve the inhomogeneous equation, one would need to 
repeat the calculations shown in the previous section, but 
this time for the whole Poincare group. To recalculate ev- 
erything with 4 more generators involves a lot of work. 
The paper by Wei & Norman ( 1963 ) provides us with a 
way to avoid some of this work. The Poincare group can be 
decomposed into the direct sum of a semi-simple algebra 
L and a radical R (whose definitions can be found in that 
same paper for instance). In terms of the previously used 
generators of the Poincare group, the semi-simple algebra 
is given by 

L = {Hi, H2, H3, Li, L2, L3}, 

the generators of the homogeneous part. The radical is 
given by 

i? = {]l, Pi,Pq,Pu,Pv}, 

the inhomogeneous part plus the identity. We will include 
the dilatation transformation in the semi-simple algebra 
set for easiness. If we write the transfer equation as 

then, H, an element of the Poincare group plus dilatations, 
can be decomposed into 

H = -K + P, 

where K is the usual absorption matrix, an element of the 
homogeneous group, and P, which stands for the set of four 
translations introduced in Section 3, is an element of both 
the radical and the inhomogeneous part of the equation. In 
the last section we dealt with the homogeneous equation 
and found a solution for the evolution operator 0(r, tq) by 
using the Lorentz group plus dilatations. Now, it is easy 
to demonstrate that if we are able to solve the equation 



a solution for the complete transfer equation can be writ- 
ten in the form 



/(T)=0(T,ro)Jfi(T). 



(26) 




'0-Ht,ro)Pit)0{t,To))lt 



In fact, this result is exactly equivalent to the for- 
mal solution given by Landi Degl'Innocenti & Landi 
Degrinnocenti(1985). To prove it, we note that 0{t, tq)I n 
will give, by properties of the evolution operator, a new 
InXt). The effect of P is however independent of the actual 
value of Ifj{t), we will always obtain that 



PI Bit) 



where the ji 's are the infinitesimal parameters of the trans- 
lation transformation: the emission vector in our particu- 
lar case. The previous equation results therefore in 

which can be integrated at once: 

Ir{t) = Ir{to) + f 0-\t, To)J{t)dt. 

J To 

Combining it with the homogeneous solution, we obtain 
the final complete solution: 

J(t) - 0(r,To) (^Ir{to)+ J\-Ht,To)J{t)d?j . 

And benefiting from the well known properties of the evo- 
lution operator, we can transform this expression into the 
formal solution given in the above referred paper: 

J(t) - 0(r,To)/fl(To) + f 0{T,t)J{t)dt. 

J To 

Hence, once the evolution operator is solved as shown 
in the previous section we can use this expression to ob- 
tain the complete solution. Instead of doing that, we shall 
proceed with the techniques provided by the group the- 
ory and obtain a completely equivalent but independent 
expression for Ir . 

P belongs to the radical which, by definition, is an ideal 
of the Poincare group, so that in fact the term 0~^P0 is 
just a linear combination of the infinitesimal generators of 
R: 

(0-ipO) = 77/ 11 + DoPi + DiPq + D2PU + D^iPv 

Obtaining the coefficients Di ( with i = 0, 1, 2, 3) is quite 
long, the detailed calculation is to be found in the Ap- 
pendix. This calculation constitutes by itself a demon- 
stration of the first statement of this paragraph for our 
particular long one, but which does not require 
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further knowledge in group theory. The next step is to 
solve the equation 

■^Ir = [DoPi + DiPq + D2PU + DsPv) Ir- 

This is in fact a very easy equation, as every commutes 
with each other. The solution can be given at once as 

Pi r D„{t)dt Pq r Di(t)dt 

I r\t) — e ° • e •^^0 

Pu r D2{t)dt Pv r D3(t)dt 

e ■>n) -e •'^0 Ir[to)- 

Calculation of the exponentials is straightforward: The Pi 
are the infinitesimal generators of translations in the four 
axes I, Q, U, V, hence by exponentiation we recuperate the 
finite transformation: 



Ir{t) = Ir{to)+ 



7fl(To)+D(r,ro).(27) 



The last step is to put together the homogeneous and 
inhomogeneous solutions by using expression ( p6|) . We ob- 
tain 

/(r) =e"^^>''*^'* []1 + 53(t)H3] • [ll + gaWHs] • 
[cosgi(T)]l + sin.gi(T)Hi] • [1 + 56(^)13] • [1 + ,g5(T)L2] ■ 

[cos54(r)]l + sing4(T)Li] • (/(ro) + D{t, to)I28) 

Note that this solution is general for all the radiative trans- 
fer problems known to date in polarization, provided the 
source function is given. As discussed in the previous sec- 
tion, it is even too general. 

This solution is also independent of any model at- 
mosphere. This is necessary to ensure its generality, but 
presents the problem of the integrability: have the inte- 
grals for the gi's and the inhomogeneous vector D an an- 
alytical expression for all and every interesting case? The 
most likely answer is no. But whatever the answer, physi- 
cal intuition indicates that there must always exist at least 
a numerical solution to them. However further work must 
be developed on the subject. 

It is also important to note that we are proposing not 
just an expression as solution of the RTE, but a method: 
particular cases may ask for different orderings of the gen- 
erators or even a different decomposition of the product of 
exponentials. We have seen that we can in any case give 
a solution in the form of seven exponentials, but, for in- 
stance, when solving constant matrix atmospheres it may 
be more interesting to consider only the product of 3 ex- 
ponentials, one with the generators a second one for 
the Li's, and a last for the dilatations, or even a sole one, 
in which case the solution for the evolution operator can 
be written at once: 

0(t,to) = exp-K(r - tq), 



in accordance with Magnus' solution (Magnus, 1954) or 
with the scalar-like exponential solution (Semel & Lopez 
Ariste, 1999). For any number of exponentials, the method 
will work, the sole problem being to solve the subsequent 
scalar linear equations and integrals. In all the cases, a 
compact and finite expression for the solution is obtained 
and the problem is reduced to the ability to integrate 
scalar expressions. 

6. Discussion and conclusion 

In this paper we have introduced a new formalism to han- 
dle Stokes parameters and radiative transfer equations for 
polarized light. In this formalism, the Stokes parameters 
appear as a 4-vector in a Minkowski-like 4-dimensional 
space, and its evolution in time looks mathematically as 
typical rotations, contractions and translations in this 
space. These movements are completely described by the 
transformations of the group of Poincare plus dilatations, 
a 10 -f 1 dimension group, well-known from other areas 
of physics and mathematics. The RTE is shown to be an 
infinitesimal transformation of this group. We therefore 
propose that a solution to the RTE can be given in the 
form of a finite transformation of the Poincare plus di- 
latations group. Obtaining of this solution from the vari- 
ables present in the transfer equation raises some techni- 
cal difficulties which have been overcome by the use of the 
Wei- Norman method (Wei and Norman, 1963). The final 
obstacle is reduced to a scalar Riccati equation. 

The Riccati equation is a well studied first order differ- 
ential equation, characterized by its quadratic term. This 
non-linearity can at worst prevent an explicit solution, 
and usually make it difficult to calculate. In any case the 
problem of giving a solution for the RTE will have been 
reduced from solving a 4-dimensional vector equation to 
solving a scalar Riccati one. Whenever this Riccati equa- 
tion can be integrated, a complete solution is obtained for 
the RTE. 

Until now only numerical integration methods (see for 
instance Rees et. al, 1989, Bellot Rubio et. al, 1998 or 
Lopez Ariste & Semel, 1999) were capable of integrating 
non-constant K matrices. The only way to test the va- 
lidity of the solution and the convergence rates was to 
compare them with previous methods, known to converge 
asymptotically. The solution presented in this paper may 
allow a comparison with an analytically exact solution. 
We anticipate that new numerical methods will be devel- 
oped taking advantage of the analytical solution; perhaps 
faster and more precise than previous ones. 

In order to obtain this solution we made use of a math- 
ematical frame, group theory, rarely seen in the astrophys- 
ical literature. The advantages gained in the integration 
of the polarized RTE warranted the efforts. We anticipate 
that new results in the study of polarized light transfer in 
astrophysical problems will be achieved by the use of this 
and related techniques. 
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Appendix A: Solution for the radical for the dilatation operator D (expressed in matrix repre- 

. ,.,.„. n , , , . , sentation by the identity matrix), we have: 

As explained m Section |5|, we need to calculate, m order 

to obtain the final solution, a term of the form jp. Q] — p. (A 2) 

(A.l) ^ith^ = /,Q,Ai?. 

where is the evolution operator, solution of the homo- Once we have all the rules of the game we can begin 
geneous equation i e P^^y with expression ( AJ ) and calculate its first term: 



0(t,to) = e- 



,ff3(r)H3 g2(T)H2 ffi(T)Hi g6(T)L3 ff5(-r)L2 g4(r)Li ff7(-r)l 



and where P is the inhomogeneous part of the transfer (in what follows and for the sake of clarity we leave out 

the dependences on r of g.; and ji to recuperate them in 
the final expressions). To this end we will need to calculate 



equation which can be written as 

P = Ji{t)Pi + jq{t)Pq + ju{t}Pu + JviT)Pv, 



where the ji are the components of the emission vector. 
In a further effort to simplify the calculations, instead of 
this linear combination we will use 

P = J I P/ + Pq + JA Pa + Jb Pb , 

where Pa = Pu + "iPv, and Ps is its complex conjugate. 
Consequently, Ja and Jb are given by 



JA 



Jb = -^Uu+ijv) 



Working out expression (A.l) implies the use of commu- 
tators of Hi, Li and the dilatations with P^. Those com- 
mutators, which can be found in any textbook on group 
theory, are, for Pj 

[P/,Hi] = -iPg 
[P/,H2] = -Pb 
[P/,H3] - -Pa. 

For Pq we have 

[Pq,Hi] = ~tPi 
[Pq,H2] = Pa 
[Pq,H3] = -Pa- 

For Pa 

[Pa, Hi] ^ iPA 
[Pa,H2] = -2Pi 
[Pa,H3] - 0. 



And for Pb 

[Pb,Hi] = 
[Pb,H2] = 
[Pb,H3] - 



-iPB 
-2Pq 

-2(P/-Pq). 



and add afterwards all the terms of the form 

Each one of which is to be calculated using an equivalent 
of the Baker-Hausdorff formula, which, for example for 
P/, affirms that 

g-93H3 p^gSsHa ^ p^ + [p^^ Hg] + ^53 [[Pi: Ha] , H3] + . . . 
The result of these calculations is 

e-93H3p^g<,3H3 ^ pj_g^p^^ 
e-93H3p^g93H3 = Pg-ggP^, 

g-S3H3p^gg3H3 ^ p^^ 

e-93H3p^gg3H3 = + 253(P/ - Pq). 
So one obtains 

Ri = {ji + 253js)P/ + Uq - 2g3jB)PQ + 
{-jigs - iQ.93 + Ja)Pa + JbPb = 
= cioP/ + ciiPg + C12PA + ciaPs- (A.3) 

The meaning of the coefficients cu is self-evident. Next 
term is 

Partial results involved are 

g-92H2p^gS2H2 ^SoPl + SiPq + S2PA + S3PB 

e-92H2p^g92H2 =SoPq- SiPi -S2PA + S^Pb 

g-S2H2p^g92H2 ^ ^^p^ ^ ^^p^ _^ 252P/ - 253PQ 
g-92H2p^g<,2H2 ^ ^.^p^ ^ g^p^ ^ 2^2Pq - 253P7, 



where the Sq, Si, S2, S3 are shortcuts for 

50 — cosh 32 • cos ^2, 

51 = sinh52 • sin 52, 

52 = (cosh 32 • sin 52 + cos 32 • sinh52), 



Commutators for the Li can be obtained as the complex 
conjugated of the corresponding ones for the Hi. Finally, S3 = (cosh 52 • sing2 - cos 52 • sinhg2). 



(A.4) 
(A.5) 

(A.6) 
(A.7) 



A. Lopez Ariste & M. Semel: Analytical solution to RTE 



11 



The result for R2 is 

R2 = (Cio^o - CllSl + 201252 - 2ci3S'3)P/ + 

(cio5i + cii5o - 2C12S3 + 2ci352)Pq 

+ (cioS'2 + CllS's + 01250 + Ci3Si)Pa + 

(cioft - C11S2 + C12S1 + Ci3So)Pb 

= C20P/ + C2iPq + C22PA + C23PB. (A.8) 
Next term is 

and, by means of the following partial results: 
g-SiHip^gSiHi ^ cOSfifiPj -isinfifiPg 
g-siHip^gSiHi ^ cosflfiPg -isinfifiP/ 

e-si^ip^efi"! = e'^iPA 
c-s^^'Pee^'^' = e-'siPs 
one gets 

R3 = (^20 cos .91 - ic2i sin.gi)P/ + 

+ (c2i cosgi - ZC20 sin 511) Pg + 

+C22e'»^PA + C23e-'»^PB = 

= C30P7 + C3iPq + C32PA + C33PS. (A.9) 

Now the process is to be repeated for to obtain R4, 
and i?6- Being the Lj the complex conjugated of Hj, every 
expression is immediate just by using the corresponding 
complex conjugated coefficients and by substituting the 
functions 54, and ge for gi, g2 and gs respectively. We 
successively obtain 

-R4 = (C30 + 2g6C33)P/ + (C3I - 256C33)Pq + 
+ (-£3056 - C3lfif6 + C32)Pa + C33PB 

= C40P7 + C4iPq + C42PA + C43PS, (A. 10) 

and 

R5 = = {c4oTo - C41T1 + 2C42T2 - 2C43T3)P/ + 

+ (c4ori + C41T0 - 204213 + 2c43T'2)Pq + 
+ (C4072 + C41T3 + C42T0 + C43ri)PA + 
+ (C40T3 - C41T2 + C42T'i + C43ro)PB 

= C50P7 + csiPq + C52PA + C53PS, (A.ll) 

where 

To = coshfii5 ■ cosfifs, (A.12) 

Ti = sinhg'5 ■ sin^is, (A. 13) 

T2 = (cosh 55 • sin 55 + cos 55 • sinhgs), (A. 14) 

^3 = ^(cosh.g5 • sin 55 - COS55 • sinh^s). (A. 15) 

The final result is 

-Re = (c5o cos .g4 + ic5i sin54)P/ + 

+ (C51 cos 34 + iC50 sin54)PQ + C52e"'^''PA + 
+C53e'»*PB = 
= C60P7 + C6iPq + C62PA + C63Ps. (A. 16) 



And we are only left with the dilatation operator, for 
which the operations are at this point almost immediate 
and give: 

(0-ipO) = e»^ (C60P7 + C6iPq + C62PA + C63Pb) • (A.17) 

The Di coefficients at section 5, can straightforwardly be 
obtained from this expression as 

Do{t) = eS^(^)c6o(T) 
Di(r) =e3^(^)c6i(T) 

£)2(r)=e9^(^)i(c62(r)+C63(r)) 

Dsir) = e»^(-)ii(c63(T) - C62(t)). (A.18) 

Acknowledgements. The authors are indebted to M. Landolfi 
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